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ABSTRACT 

The large scatters of luminosity relations of gamma-ray bursts (GRBs) have been one of the 
most important reasons that prevent the extensive applications of GRBs in cosmology. In this 
paper, we extend the two-dimensional (2D) luminosity relations with rjag, V, -Epcak: and trt as 
the luminosity indicators to three dimensions (3D) using the same set of luminosity indicators to 
explore the possibility of decreasing the intrinsic scatters. We find that, for the 3D luminosity 
relations between the luminosity and an energy scale (-Epoak) and a time scale (rjag or trt), their 
intrinsic scatters are considerably smaller than those of corresponding 2D luminosity relations. 
Enlightened by the result and the definition of the luminosity (energy released in units of time), 
we discussed possible reasons behind, which may give us helpful suggestions on seeking more 
precise luminosity relations for GRBs in the future. 

Subject headings: Gamma rays: bursts 



1. Introduction 

Gamma-ray bursts (GRBs) have recently at- 
tracted much attention in their cosmological ap- 
plications as the most luminous astrophysical 
events observed today. Based on the correla- 
tions between the luminosity/energy and the 
measurable parameters of light curves and/or 
spectra, GRBs can be used as stan dard can- 
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The advantage of the GRBs is their high red- 
shifts due to the ir high lumino sities. The 69 



GRBs compiled in lSchaefen (|2007f ) extend the red- 
shift to z > 6. Recently obs erved GRB 090423 
has a redshift of z = 8.3 ( Tanvir et al. 120091: 
Salvaterra et al.ll2009h . However, GRBs are not 



as ideal standard candles as type la supernovae 
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(SNe la). Compared with SNe la, GRBs suffer 
from the circularity problem due to the lack of 
low-redshift samples and the scatters of known 
luminosity relations of GRBs are still very large. 
There are a few ways proposed in literatures to 
avoid the circularity problem. One can simply fit 
the calibration paramete rs and cosmological pa- 
rame ters simu ltaneously (iLi et al. 2008t lOi et al 
20081 ): and in IWand (l2008|), GRB data are sum- 



marized by a set of model-independent distance 
measurements; calibrating GRBs using SNe la 
in thei r overlapping redshift range was also pro- 
posed (iKodama e t al. 2008; Li ang et al. 12008') and 
adopted (|Wei k Zhana.200a : .Cardone et al.ii2009l ) 
in cosmological studies using GRBs. 

Till now in most works of cosmological stud- 
ies using GRBs, two-dimensional (2D) luminosity 
relations are used, which have leas t calibration 
param eters. For example Tia.s-L ( Norris et al.l 
2000l). V-L ("Feni more fc Ramirez-Ruid l200d 
Reichart et a l. 2001, there exist several defini- 



tions of V, mainly depending on the smooth- 
ing time intervals the reference curve is built 
upon, and on the norma lizat:' " 
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20071) relations. However, the intrinsic scatters of 
2D luminosity relations are usually very large, 
which may imply hidden parameters considering 
the complication of GRBs. There are already 
works which explored the possibility of a three- 
dimensional (3D) correlation with negligible scat- 
ter. For example, iFirmani et al.l (j2006t ) claimed 
that a temporal parameter of the prompt emis- 
sion, the To. 45, could reduce the scatter of the cor- 
relation of LisQ-Epcak to a negligible value. But 
it was later found that the new proposed relation 
doe s not appear to be as tight as it seemed t o 
be (jRossi et all 120081 : ICoUazzi fc Schaefed liooi ) . 
In this paper, we extend the luminosity relations 
used in Schacfcr (2003) from 2D to 3D and ex- 
plore the possibility of decreasing the scatters in 
the correlations. 

2. Methodology 



In ISchaefeij (|2007r ). five luminosity relations of 
GRBs as follows are used: 
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where the luminosity L and the total coUimation- 
corrected energy of GRBs are derived respec- 
tively from the bolometric peak flux Pboio and the 
bolometric fluence Sho\o of GRBs through 



(6) 



E^ — i^^^so-Fbcam — 47rd|S'boio(l + ^FbclZ) 

Here di is the luminosity distance, which depends 
on the cosmological model and is inversely propor- 
tional to the value of Hubble parameter of today. 
In this paper, we have adopted the flat ACDM 
model with = 0.27 and in the calculation, we 
actually replace with J/, = -^d^ x 1 cm in 
Eq. © and Eq. ([7]), so that the dependence on 
Hubble constant is absorbed into the intercepts of 



the linear luminosity relations. For later conve- 
nience, we denote these luminosity relations by 



where 
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The coefficient c is given two superscripts to in- 
corporate 3D correlations introduced just below. 
Hereafter, we denote a luminosity relation by the 
superscript pair of the corresponding c. 

The above 2D luminosity relations connect 
measurable parameters of light curves and/or 
spectra with GRB luminosity /energy. They are 
emp irical though th ere exist some explanations 
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Since the physical processes of GRBs are very 
complicated, the real correlations between GRB 
luminosity/energy and the parameters of light 
curves and/or spectra, generally speaking, should 
be more complicated than the above simple rela- 
tions. The large scatters in the above luminosity 
relations also imply that there may be hidden 
parameters not included. Motivated by this, we 
extend the 2D luminosity relations to 3D and 
investigate if there are any improvements. See 
Eq. dlT]), Eq. ([13]), and Eq. since a;(3) = a;^ 

and y^'^^ is different from y^^^ {i — 1, 2, 3, 5), we 
introduce 3D correlations as follows. For (i, j) 
with both i and j in (1, 2, 3, 5) and i < j, we 
introduce 3D correlations as 
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The luminosity relation (3, 4) is different from 
other 3D luminosity relations in that there is only 
one luminosity indicator (i.e. i^pcak) in it, while 
there are two in other 3D luminosity relations. 
The luminosity relations in Eq. (|16p and those in 
Eq. (fT8|) are different from each other in that the 
left-hand side of Eq. p6|) is log ^ ^^g^-i , while the 

left-hand side of Eq. HID is log 

We investigate the quality of the luminosity re- 
lations mainly by comparing their intrinsic scat- 
ters. Since we are free to multiply the luminos- 
ity relations by a constant when introducing 3D 
correlations and such a factor would increase the 
intrinsic scatter by the same multiple as the fac- 
tor, we need to normalize the relations in order to 
compare the intrinsic scatters. Bearing in mind 
that one of the most important purposes for ex- 
ploring the luminosity relations is to use them in 
distance measurements, what we do here is just 
dividing Eq. pT)) by a factor of 1 — c^' so that 
log(dL) have the same coefficient in all luminosity 
relations above. 

Comparing the 3D luminosity relations with 
the 2Ds, as an example, comparing Eq. p6p with 
Eq. (HI), one can see that, by introducing the 3D 
correlations in this way, we actually treat x^^^ 
as the hidden parameter of the 2D correlation of 
Eq. ^ and write it explicitly in the 3D correlation 
of Eq. (|16p . So, in addition to examine the qual- 
ity of the luminosity relations by comparing their 
intrinsic scatters, we also calculated the correla- 
tion coefficients between the residual of the fit of 
2D correlations and the possible hidden parame- 
ters we introduced to extend correlations from 2D 
to 3D. 

In the fit of the lumino sity relations, we used 
the techniques presented in lD'Agostinil (j2005[ ). fol- 
lowing which the likelihood function for the coef- 
ficients c and the intrinsic scatter is (for the cases 
of Eq. (|16p . The other cases of 3D luminosity re- 
lations are similar and the likelihood function for 
the 2D luminosity relations can be obtained just 
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where k runs over GRBs with corresponding quan- 
tities available. In the calculation, Markov chain 
Monte Carlo techniques are used. For each lumi- 
nosity relations, A Markov chain with samples of 
order 10^ is generated according to the likelihood 
function and then properly burned in and thinned 
to derive statistics of interested parameters. 

For the GRB d ata, we used the compilation 
in ISchaefeil (|2007h . which includes 69 GRBs. 



When considering error propagation from a quan- 
tity, say ^ with error cr^, to its logarithm, we set 
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center value and the error of the logarithm corre- 
spondingly. This requires ^ > (the quantities 
we are interested in here are all positive). Due to 
the limitation of the data, for a given luminosity 
relation (i, j), not all the GRBs have all of the 
needed observational quantities available and sat- 
isfy ^ > (T^ at the same time. By set («, j) we 
denote the maximum GRB set that can be used 
in the luminosity relation (j, j). The numbers of 
GRBs of different sets are presented in Table [1] 

3. Results and discussion 

We summarize our results in Tables [T] and [2j 
In Table [U by comparing the intrinsic scatters of 
the 3D luminosity relations with those of the cor- 
responding 2D luminosity relations, we find that 
only for the cases of (1, 3) and (3, 5) the intrinsic 
scatters of the 3D luminosity relations are consid- 
erably smaller than those of their corresponding 
2D luminosity relations, and for all other cases, 
the intrinsic scatters of the 3D luminosity relations 
are either very close to the smaller one of the in- 
trinsic scatters of the corresponding 2D luminosity 
relations (for correlations in Eq. (fT8|) , close to the 
intrinsic scatter of luminosity relation (4, 4)) or 
even greater than those of the corresponding 2D 
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Table 1: Fit of 2D and 3D luminosity relations. In every grid, the first row is the number of GRBs of set (i, j), the vector below enclosed 
by parentheses is the vector of c for the luminosity relation (z, j), and what follows next is the intrinsic scatter. For 3D luminosity relations, 
their reduction in the intrinsic scatters compared to corresponding 2D luminosity relations are presented in the brackets: for 3D luminosity 
relations in Eq. ([TBI) and Eq. ([17]), the reduction corresponds to the 2D luminosity relation (z, i) and (j, j) in turn and for those in Eq. (jlSp. 
corresponds to (4, 4). The statistics in the table are for the median values and the errors of la (68.3%) confidence level. 
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Table 2: Correlation coefhcients between the residuals of the fit of 2D correlations and possible hidden parameters. The statistics in the 
table are for the median values and the errors of la (68.3%) confidence level. 



luminosity relations (for the case (3, 4)). The cor- 
relations presented in Table [2] also give the same 
implication. Only the correlation coefficients be- 
tween the residual of fit of (1, 1) and x^^\ (3, 3) 
and x^^\ (5, 5) and a;'-^-' are approximately greater 
than 0.7. We presented the plots of the luminosity 
relation (1, 3) and (3, 5) in Figure[TJ 





=1 log -^2 106=^ 



Fig. 1. — Plots of the 3D luminosity relation (1, 3) 
and (3, 5). 

It is interesting to note that the 3D luminosity 
relations that have considerable improvement in 
the intrinsic scatter are the luminosity relations 
in Eq. pop . with one of the luminosity indica- 
tors being an energy scale (i?pcak) and the other a 
time scale (riag or trt)- Enlightened by this, one 
may naturally guess that probably the luminosi- 
ties of GRBs are mainly determined by a charac- 
teristic energy scale and a characteristic time scale 
and i?poak is correlated to the characteristic energy 
scale, Tiag and trt to the characteristic time scale, 
so that 3D luminosity relations with Epcak and riag 
or Trt as luminosity indicators could significantly 
reduce the intrinsic scatters compared to corre- 



sponding 2D luminosity relations. This is easy to 
understand, since the definition of the luminos- 
ity is the amount of the energy released in units 
of time. If one knows the released energy and the 
time duration of a GRB, the (averaged) luminosity 
can be calculated immediately. Obviously, the 2D 
correlation between the luminosity and the energy 
or the time duration would not be complete unless 
the time duration or the energy is constant for all 
samples or there is a strong correlation between 
the energy and the time duration. This, if proven 
true, also explains why the 2D luminosity relations 
(3, 3) and (4, 4) have little correlation with each 
other despite their simil arity at fir s t sigh t (see dis- 
cussion in Section 4.7 of Schaefer (|2007l) : the cor- 
relation coefficient between the residual of fit of 
(3, 3) and y^''^ and the correlation coefficient be- 
tween the residual of fit of (4, 4) and y^^^ also show 
their weak correlation). This is because of the 
difference between and L, where a new inde- 
pendent variable — the characteristic time scale — 
enters. In fact, we could check the guess partially 
with current data by examining the correlation be- 
tween a;(*)s. InTableEl one can find that cc^^^ , x'^-* , 
and x^^^ are strongly correlated with each other, 
while all of th em have only w e ak co rrelation with 
x^^') ( see also ISchaefer et al. ( 2001 ) and Section 



4.5 of ISchaefej (|2007h for discussions on the cor- 



relations between the luminosity indicators) . This 
is consistent with the guess that the luminosity 
indicators Tiag and trt are correlated with a char- 
acteristic time scale and -Epcak is correlated with a 
characteristic energy scale which is independent of 
the characteristic time scale. In addition, the cor- 
relations presented in Table [3] seems to imply that 
V is also correlated with the characteristic time 
scale. Accordingly, there is indeed some reduction 
in the intrinsic scatter for the 3D luminosity re- 
lation (2, 3) compared with corresponding 2D lu- 
minosity relations, but such reduction is relatively 
smaller than that of (1, 3) and (3, 5). Maybe the 
correlation of V with the characteristic time scale 
is not so strong as that of Tjag or trt- 

If our guess about the GRB luminosity re- 
lations is correct, it would be very enlighten- 
ing. It suggests us to include an energy scale 
and a corresponding time scale for seeking more 
precise luminosity relations for GRBs in the fu- 
ture. However, it should be emphasized that, even 
if it is true, only appropriate energy and time 
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scales might significantly reduce the intrinsic scat- 
ter bearing in mind the situation of some 3D lu- 
minosity relations (see, for example, Firmani et alj 
2006HRossi et allbOOst ICoUazzi fc Schaefer..2008lV 



4. Summary 

In this paper, we extend the widely used 2D 
luminosity relations to 3D by using the same set 
of luminosity indicators, i.e. riag, V, iJpoak, and 
trti and check the improvement in the quality of 
the luminosity relations. We find that, for the 
3D luminosity relations between the luminosity 
and an energy scale (iJpeak) and a time scale (riag 
or trt), their intrinsic scatters are considerably 
smaller than those of corresponding 2D luminosity 
relations. The correlations between the residuals 
of fit of the 2D luminosity relations and the lumi- 
nosity indicators also give the same implication. 
Enlightened by the result and the definition of the 
luminosity (energy released in units of time), we 
discussed possible reasons behind, which may give 
us helpful suggestions on seeking more precise lu- 
minosity relations for GRBs in the future. 
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